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Abstract. The present article deals with differential equations 
with spectral parameter from the point of view of formal power se- 
ries. The treatment does not make use of the notion of eigenvalue, 
but introduces a new idea: the spectral residue. 

The article focuses on second-order, self-adjoint problems. In 
such a setting every potential function determines a sequence of 
spectral residues. This correspondence is invertible, and gives rise 
to a combinatorial inversion formula. Other interesting combinato- 
rial consequences are obtained by considering spectral residues of 
exactly-solvable potentials of 1-dimensional quantum mechanics. 

It is also shown that the Darboux transformation of 1-dimensional 
potentials corresponds to a simple negation of the corresponding 
spectral residues. This fact leads to another combinatorial inver- 
sion formula. 

Finally, there is a brief discussion of applications. The topics 
considered are enumeration problems and integrablc systems. 
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1. Introduction 



There is a rich interplay between the methods and ideas of theoreti- 
cal physics and that of combinatorics. The connection runs both ways: 
combinatorial methods are applied to the study of physical models, 
while ideas from physics are used to obtain results of a purely combi- 
natorial nature. 

The present paper follows the second pattern. The aim is to demon- 
strate that certain methods for the solution of the time-independent 
Scrhodinger equation lead naturally to interesting identities and inver- 
sion formulas. The idea that allows us to go from the one domain to 
the other is that of a spectral residue. The application of differential 
operators and residue techniques to combinatorial identities and inver- 
sion formulas is well documented [l], 0, 0, (3). However, the notion of 
a spectral residue seems to be unknown, and therefore deserves to be 
considered in some detail. 

The impetus for spectral residues comes from the following problem: 
to give a formal series treatment of differential equations with spectral 
parameter. It is far from clear whether the notion of eigenvalue is 
meaningful in a formal setting, and so we take a different approach. 
Consider a linear differential operator T in one independent variable 
z, with spectral parameter v, and define 4>(z; v) as the formal series 
solution of T[4>] = 0. As the analogue of the spectrum of T, let us take 
the set of values of v for which some of the coefficients of are singular. 
This spectral set is not of interest in itself; in the cases considered here 
it is just the set of natural numbers. More interesting are the residue 
of the coefficients of the solution at the singular values of v\ these we 
will call the spectral residues of the operator. 

Let us restrict our attention to classical quantum-mechanics and con- 
sider second-order, differential equations that are amenable to the usual 
method of undetermined coefficients (aka the method of Frobenius): 

T[4>\ = z 2 <p"{z) + (1 - v + P{z))z<j)'{z) + Q{z)<P{z) = 0, 0(0) = 1. 

Here P and Q are formal power series in z with vanishing constant term, 
and v is the spectral parameter. The roots of the indicial equation are 
and v. Hence if v is a natural number, then, in general, the solution is 
singular. In this way there arises a sequence of spectral residues, one for 
each natural number. If we can solve the operator equation T[0] = 0, 
then we will be able to obtain information about this sequence. 

From the very beginning of quantum mechanics, the study of Schro- 
dinger's equation was marked by the search for exact solutions. The 
two outstanding methods in this endeavor are first, the theory of hyper- 
geometric functions 0, and second the Darboux transformation || [7|]. 

2 



SPECTRAL RESIDUES 



3 



The first method furnishes a large class of exactly-solvable potentials 
whose eigenfunctions can be given in terms of either the hypergeo- 
metric or the confluent hypergeometric functions^]. All the classical 
exactly-solvable potentials: the Harmonic oscillator, the Coulomb, the 
Morse, the Eckart potentials, are of this type. 

As for the second method, the essential feature of the Darboux trans- 
formation is that it changes one potential function into another in such 
a way that the eigenfunctions of the one can be explicitly constructed 
from the eigenfunctions of the other. By means of this pairing exactly 
solvable potentials yield new exactly solvable potentials. Combining 
the two methods it becomes possible to describe a very large class of 
exactly solvable potentials. 

Our goal is to examine these two methodologies from a combinatorial 
standpoint by using spectral residues. In this context the following 
questions are quite natural. 

1. Does the spectral residue sequence completely determine the po- 
tential? If yes, can this relationship be described explicitly? 

2. What are the spectral residue sequences corresponding to poten- 
tials of hypergeometric type? Does the exact-solvability of these 
potentials manifest itself at the level of spectral residues? 

3. What is the effect of the Darboux transformation on the corre- 
sponding spectral residue sequence? 

We will see that all of these questions have satisfactory answers, and 
that these answers manifest themselves as a novel class of summation 
identities and inversion formulas. The sums that occur in this context 
are indexed by compositions, i.e. ordered partitions. More specifically, 
we will say that an /-part composition of a natural number n is an 
ordered list of I positive integers p = (pi,P2, ■ ■ ■ ,Pi) such that 

Pi + . . . + pi = n. 

On occasion we will write |p| instead of n, and £(p) instead of I. 
Identities that feature sums indexed by partitions are not new || [T(], 



e.g. §60]||11||. All of these references feature sums whose terms are 
symmetric functions of the partition entries. Thus, no generality is 
lost by working with partitions rather than compositions. By contrast, 
the sums that arise in connection with spectral residues have terms 
that are functions of the partial sums of the composition entries; the 
order of the entries cannot be ignored. 

For a composition p = . . . ,pi) we define 



Sp — S1S2 • • • Si, 
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where Sj — p\ + . . . + Pj is the j th left partial sum corresponding to p. 
Thus, if one interprets p as a walk from to n with steps Pi, ■ ■ ■ ,pi, 
then s p is the product of numbers visited on that walk. Let p' denote 
the reversed composition (pi, . . . ,p 2 ,pi). In this way s p / is equal to 
the product of visited numbers when one interprets p as a reverse walk 
from n down to 0. 

Let U\,U2,Uz,... be an infinite sequence of non-commuting vari- 
ables. For n = 1, 2, . . . set 

Pn=J2 7 ^-U p (1.1) 

i i P P 
|p[=n 

where the index runs over all compositions p = (pi, ... ,p{) of n, and 
U p is an abbreviation for the /-fold product U Pl U P2 . . . U PV Here are the 
first few elements of the resulting sequence: 

Pi = Ux, 

92 = \u 2 + l -U,\ 

p 3 = \u* + \{U 2 U X + UxU 2 ) + ^-U^, 
6 b 12 

p 4 = 7C/4 + ^(Wi + C/1C/3) + ^u 2 2 + 

4 12 lo 

+ ^(Ui 2 U 2 + U 2 U X 2 ) + ^umk + 

Evidently, equation (|1 . 1| ) defines an invertible, nonlinear mapping from 
the set of sequences {U n } to the set of sequences {p n }. 

Theorem 1.1. The inverse mapping is given by 

11 (Pi + P2) (P2 + Pa) • • • (Pi-i + Pi) 

\p\=n 

where the sum is taken over all compositions p = (pi, ... ,pi) of n, and 
Pp = PpiPp2 ■ ■ ■ Ppr 

What does this theorem have to do with potential functions and 
spectral residues? Letting U(z) = U\Z + U 2 z 2 + . . . and working with 
the "nearly" self-adjoint equation 

*V' (*) + (1 - v)ztf{£) + U{z)i){z) = 0, 

we will see below that the corresponding spectral residues are given 
by equation ( |1 . 1| ) . The inversion formula given in Theorem |1.1| is mo- 



tivated by question 1 above. The proof of the formula will be given 
below in Section 0. 
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Before continuing, we should remark that in classical quantum me- 
chanics the potential function U(z) is scalar-valued, and hence the 
coefficients U n commute. Likewise the usual application of the method 
of Frobenius is to series with commuting coefficients. However, formula 
such as and ( |1.2| ) involve a sum indexed by compositions, and as 
such are best treated in a non-commutative setting. Indeed, if the U n 
commuted, then the natural form for the inversion formula would be 
a sum indexed by partitions, i.e. unordered lists of natural numbers, 
rather than compositions. 

The reader should therefore be warned that our discussion of spectral 
residues will involve both the commutative and the non-commutative 
settings. In as much as commutativity is an extra assumption, we 
will try to do as much as possible with non-commuting coefficients. 
However, much of the discussion, in particular question 2 requires the 
commutativity assumption. As a consequence, we will find ourselves 
jumping several times between the two different settings. 

As regards question 2, we should first recall that the hypergeometric 
equation has 3 parameters. Making a suitable change of the parameter 
variables one of them can be made into a spectral parameter. Taking 
spectral residues "uses up" this parameter, and as a result the spectral 
residue sequence consists of polynomials in the remaining two param- 
eters. 

We should also observe that the exact-solvability of the hyperge- 
ometric potentials is really due to the fact that the hypergeometric 
equation admits polynomial solutions for certain discrete values of its 
parameters. This observation has an important consequence: it turns 
out that the spectral residues of hypergeometric potentials factorize 
with respect to the two remaining parameters. This in turn leads to 
interesting combinatorial identities. 

Here is an example of this phenomenon associated to the potential 
function 

U ( z ) = m — = vz + 2vz 2 + 3vz 3 + ... , 
(1-z) 2 

where the v is a parameter. (N.B. The other hypergeometric parameter 
was dropped to simplify this example. It will reappear in Section |3|, 
where this sort of factorization is discussed fully.) As per the 
spectral residues are given by 

7T7 (n*H (L3) 

tip tip' \ / 

|p|=n 
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where the sum is taken over all compositions p = (pi, . . . ,pi) of n. 
Here are the first few of these polynomials: 

Pi=v 

1 2 
P2 = V + -V 

2 1 

P3 = V + -V +-V 

Pl = v + -v +^v> + m v 

Routine calculation indicates that the above polynomials can be fac- 
tored. In fact, we have the following general result. 

Theorem 1.2. The polynomials p n (v) factorize completely over the 
rationals. Indeed, 

n\(n - 1)! p n = v(v + 1 • 2)(v + 2 ■ 3) . . . (v + (n - 1) ■ n). (1.4) 

This result can also be appreciated in its converse formulation: the 
coefficient of v k of the Pochammer-symbol-like polynomial in the right- 
hand side of (O) is given by the composition sum 



\p\=n ~f °P i=l 

e(p)=k 



where now the index ranges over only the /c-part compositions of n. The 
proof will be given in Section |[ For now let us stress that the reason 
underlying this factorization phenomenon is the exact solvability of the 
corresponding potential U(z). We should also note that this particular 
U (z) is just a special case of the well-known Eckart potential. 

As regards question 3, for now let us limit ourselves to saying that 
at the level of power series the Darboux transformation is given by a 
certain invertible, non-linear mapping from a power series 

U(z) = U x z + U 2 z 2 + ... 

to a power series 

U(z) = U lZ + U 2 z 2 + . . . . 
The mapping in question is described by a certain composition sum 

Un = W (P)^P' 
|p|=n 

where the weight w(p) is a rational function of the composition ele- 
ments Pi, ■ ■ ■ ,pi that is too complicated to be written down explicitly. 
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What is interesting, however, is how easily the Darboux transformation 
can be described in terms of spectral residues. 

Theorem 1.3. LetU andU be formal powers series with zero constant 
term. If U and U are related by a Darboux transformation, then the 
corresponding residue sequence of one is the negative of the other. 

The proof will be given in Section |j. An interesting corollary of 
this theorem (the proof will again be deferred till Section [| ) is the 
following inversion formula. For each composition p = (px, ■ ■ ■ ,P2k+i) 
of odd length set 

Qp — S2S4SQ . . . S2k, 

where again Sj — p\ + . . . + pj. Next, let pi,P2, ■ ■ ■ be an infinite 
sequence of non-commuting indeterminates, and for each n = 1,2,... 
set 

Wn=J2 , - v - \ - ( — ^P P . (1-6) 

^ (Pi +P2){P2 +P3) ■ ■ ■ [Pl-1 +Pl ] 



p| — n 
Hp) odd 



Corollary 1.4. The inverse mapping is given by 

Ee(p)-i 1 
(-1)^—-^ (1-7) 

IpI=„ 1p' 

l(p) odd 

where W p is an abbreviation for W Pl W P2 . . . W Pl . 

The organization of the remainder of the paper will mirror the above 
introductory discussion. In Section^ we introduce spectral residues, de- 
scribe how a potential function determines a spectral residue sequence, 
and give the proof of the inversion formula contained in Theorem |1 . 1| . 
In Section |3| we will treat the spectral residue sequences of the exactly- 
solvable, hypergeometric potentials, and derive number of factorization 
results similar to Theorem [L^. These results are developed in response 
to question 2 above, and serve as evidence to the assertion that exact 
solvability of a given potential function should manifest as factorizabil- 
ity of the corresponding spectral residues. Section £| we will discuss 
the Darboux transformation from the point of view of spectral residue 
sequences. The key contents are the proof of Theorem |1.3| , and of the 
inversion formula of Corollary |1.4| . In the final section we briefly survey 
applications of these identities and inversion formulas to problems in 
enumeration and to the theory of integrable systems. 
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2. Spectral residues of Schrodinger operators 

We will begin this section by recalling the details of the classic 
method of Frobenius for power series solutions of a second-order differ- 
ential equation: 

z 2 (j)"(z) + P{z)z<P\z) + Q{z)4>{z) = 0. 

We need to generalize the classical method somewhat, and therefore 
allow P(z) and Q(z) to be power series with coefficients in an arbitrary 
algebra — not necessarily commutative — over the complex numbers. 
We want, however, to retain the notion of an indicial equation, and 
therefore demand that the constant terms of P and Q be complex 
numbers. It is then easy to show that there exists at least one formal 
solution of the form 

z v {I + <p lZ + <p 2 z 2 + ...), 

where r = v is a root of the indicial equation 

r 2 + (P - l)r + Qo = 0. 

It is well known that there may not be a second such solution if the 
solutions of the indicial equation differ by an integer. 

Since at least one formal solution always exists, no generality will 
be lost if we assume that one of the roots is zero. We will therefore 
restrict ourselves to equations 

T[0] = 0, 

where T is a differential operator of the form 

T = z 2 d zz + (1 - v)zd z + P{z)zd z + Q(z), (2.8) 

and both P(z) and Q(z) are power series with zero constant term. This 
way the indicial equation is r(r — v) = 0. 

If v is not a positive integer, then there is a unique formal series 
solution <f)(z) with the initial condition 0(0) = 1. The coefficients are 
specified by the following recursive relations: 

n— 1 

n(v - n)<p n = Q n + ^2(jP n -j + Q n -j)<Pj- 

3=1 

They can also be described in terms of a composition sum: 
^ = y Sl ^P Pl + Q x x Sl P 2 + Q x Q P1 

f-^ n(is-n) s 2 {v - s 2 ) Si(z/-si) 

|p|=n 

where the index ranges over all compositions p = . . . ,pi) of n, and 
as before Sj — p± + . . . + pj. 
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When v is a positive integer, then in general there is only one so- 
lution; this solution is a series with leading term z v . However, under 
certain circumstances it is possible for a monic series solution to exist as 
well. To describe when this happens we make the following definition. 

Definition 2.1. Let T be a second- order differential operator as given 
in fl2.8|) ; but with v denoting an indeterminate rather than a complex 
number. Let <f> n be the coefficients of the corresponding formal monic 
series solution. For n a positive integer, we define the n th spectral 
residue of T to be the residue of 4> n at v = n. 

Proposition 2.2. Let T be as above, and let n be a positive integer. 
A monic series solution of 

T[0] = 

u=n 

exists if and only if the n spectral residue of T is zero. 
Proof. The n th spectral residue vanishes if and only if 

n— 1 

Qn + y~]{jP n -j + Qn-j)<Pj = 
3=1 

If this is the case, then evidently there exists a monic series solution 
for every possible choice of <p n . □ 

Another interesting property of spectral residues is their invariance 
with respect to gauge transformations. In the present context we define 
a gauge transformation to be the conjugation of a differential operator 
by an invertible multiplication operator. In particular, given T of the 
form shown in ( J2.8|) , and a monic power series = 1 + \i\z + ii^z 1 + . . . 
we obtain another such second-order operator by setting 

T = ji o T o jj,^ 1 , 

= z 2 d zz + (1 - v)zd z + P{z)zd z + Q{z). 

Gauge transformations arise naturally in the context of linear differ- 
ential equations. Indeed, if <p(z) is a monic solution of T[(f)] = 0, then 
<fi = fi<f) is evidently a monic solution of T[(j)] = 0. 

It is also important to note that the presence of the (1 — v)zd z term 
in T will cause the coefficients of Q(z) to depend on the spectral pa- 
rameter v. In the subsequent discussion we will therefore broaden our 
class of operators T so as to include the possibility that the coefficients 
of P(z) and Q(z) depend polynomially on v. 
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Proposition 2.3. If two second- order differential operators of the form 
shown in (|2.8 ) are related by a gauge transformation, then all of their 
spectral residues are equal. 

Proof. Let T, T and /i be as above, and let <p(z) and <p(z) be the 
corresponding monic series solutions. From = fMp it follows that 

<Pn — 4>n + a linear combination of 0i, . . . n _i. 

From (|2.9|) it is evident that the residues of cf>i, . . . , 4> n -i at u = n are 
all zero. The desired conclusion follows immediately. □ 

A particularly important type of second-order differential operator 
is the class of Schrodinger operators (the term Hamiltonian operator 
will also be used.) In the present setting this will be taken to mean an 
operator of the form 

H = z 2 d zz + (l-v)zd z + U(z), (2.10) 

where the potential function U (z) = U\z + U2Z 2 + . . . is a formal series 
with vanishing constant term, and where again we allow the coefficients 
Ui, U2, ... to depend polynomially on v. 

It is more customary to define a Schrodinger operator as a formally 
self-adjoint operator of the form 

H = z 2 d z + zd z + U(z) - A, (2.11) 

where A is the spectral parameter, and z is related to the physical 
distance variable x by z = e x . However, the formal eigenfunctions of Ti 
are of the form z ±y ^ times a monic series. To get at the monic series 
directly we conjugate Ti by z v l 2 , with A = ^ 2 /4, to obtain an operator 
in the "nearly" self-adjoint form shown in ( |2.10|) . 

It is important to note that the class of Schrodinger operators pro- 
vides a canonical form for second-order differential operators with re- 
spect to gauge-transformations. 

Proposition 2.4. Let T be a second-order operator of the form shown 
in ( |2.8| ). There exists a monic series /j,(z) and a potential function U(z) 
with zero constant term such that 

z 2 d zz + (1 - v)zd z + U(z) = fioTo fT 1 . 
Proof. In the non-commutative setting it suffices to take 

if 71 ' 1 



3=1 

U = (fj,Q — z 2 fi" — (1 — ijT , 
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where P(z), Q(z) are the coefficient series of the operator T. In the 
commutative setting one can use simpler formulas. Indeed, it suffices 
to introduce an auxiliary function 

a{z) = I ^T' a(0) = ' 

and then set 

/i = e CT , 

U = Q - z 2 a" - (1 - v)zo' - z 2 (a') 2 . (2.12) 

□ 

Gauge invariance of the spectral residues will be of crucial impor- 
tance when we consider hypergeometric potentials in the next section. 
For the remainder of the present section, we will focus our attention 
on the spectral residues of a Schrodinger operator. 

Proposition 2.5. The n th spectral residue, call it p n , of a Schrodinger 
operator with potential U(z) = U\z + U 2 z 2 + . . . is given by 

n 



P n = Y] Up 



Sp Sp' 
\p\=n 



Proof. Taking P(z) = and Q(z) = U(z), the desired conclusion fol- 
lows immediately from (|2.9|) above. □ 



From formula ( |1 . 1|) we see that p n = n~ 1 U n + lower order terms. 
It is therefore evident that the sequence of spectral residues p%, p2, ■ ■ ■ 
completely determines the potential. Surprisingly, the formula for the 
inverse mapping can be given explicitly. This is formula (|1.2j) of The- 
orem |1.1| . We conclude the present section by giving a proof of this 
theorem. 

Let H be a Schrodinger operator of the form shown in (|2.10 ), and 



let <f)(z) = 1 + <f>\Z + <p2Z 2 + . • • be the formal monic series solution of 

H[<t>] = 0. 

Specializing (|2.9|) we see that 



0n=E^S X --- X ^ X ( UP1 V ( 2 - 13 ) 

f-^niy-n) s 2 [y - s 2 ) si(v - si) 

\p\=n 

For all pairs of positive integers k and n with k < n, we let p n ^ denote 
the residue of <f) n at v — k. Thus, p n ^ n is just another way to refer to 
the n th spectral residue p n . 
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We are now able to write 



Ph. a _|_ Pn,n—\ _|_ _|_ Pn,l ^ 



v — n v — n + 1 v — 1 

Next, we evaluate p ni fc in terms of the coefficients of U(z). From 
( |2.13| ) we obtain 

= Uqm ■ ■ ■ Ugi Upi - ■■ Upi 

(k + t m )(-t m )...(k + t 1 )(-t 1 ) k...(k-s 2 )s 2 (k-s 1 )s 1 , 

I P I ™ 
|q|— Ti — A: 

where the sum is indexed by all compositions p = (pi, . . . ,pi) of k and 
all compositions q = (gi, . . . , q m ) of n — k, and where Sj — p\ + . . . +pj 
and tj = qi + . . . + qj. It is evident by inspection that the above sum 
can be factored, and indeed that 



Pn,k (pn—k 

Using ( 2.14j) once again we obtain 



x p k - 

-k 



I Pn—k , Pn—k,n—k—l , , Pn—k,l \ /r> 1 r\ 

pn > k = -[kT^ + TT^T + --- + TTT) pk - (2 ' l5) 

Now is a solution of TC[(f)] =0 if and only if 

n-l 

n(v - n)(f) n = U n + ^2 u n-j(pj- 

3=1 

Substituting (|2.14j ) into the left hand side, and evaluating both sides 



at v = oo we obtain 

U n = n (p n + p n ,n-l + • • • + Pn,l) ■ 

This equation together with ( |2.15|) gives a recursive specification of U n 
in terms of p\, . . . , p n . Rewriting this as a composition sum we obtain 
the desired formula (11.21). □ 



3. Spectral residues of exactly solvable potentials 

We have already mentioned in the introduction that the spectral 
residue sequence is sensitive to the exact solvability of the underlying 
potential. In the present section we will discuss this phenomenon in 
some detail, and use it to derive a number of factorization identities 



similar to the identity described in Theorem [L2] above. 

The key to the whole matter is the invariance of spectral residues 
with respect to gauge transformations — see Proposition ^73] above. 
Hypergeometric functions also play an important role. Indeed the ex- 
act solvability of the potentials that we will be considering here is based 
on the fact that the corresponding Schrodigner operators are gauge 
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equivalent to the hypergeometric operator. Certain values of the hy- 
pergeometric parameters cause the hypergeometric series to truncate, 
which means that for those values of the parameters almost all the 
spectral residues vanish. Thus, when we consider spectral residues as 
polynomials in the hypergeometric parameters we obtain factorizations 
of said polynomials. 

Let us also remark that in the present section we can safely restrict 
ourselves to power series with complex coefficients. We will consider 
non-commutative coefficients again in section [|. 

Three families of exactly-solvable potentials will be treated here: the 
Eckart [15], the Poschl- Teller [IB], and the Morse [TJ| potentials. The 



first two of these have eigenfunctions related to Gaussian hypergeomet- 
ric functions; in the latter case one needs to use confluent hypergeo- 
metric functions. 

We will discuss each of these families in turn. In each case we will 
first show how the Schrodinger's equation is related to the hypergeo- 
metric or the confluent hypergeometric equation. Next we will consider 
the corresponding spectral residues and see that in each case there is 
complete factorization in terms of the potential parameters. Finally, 
we will list some interesting factorization identities that arise for each 
of the families for certain specialized values of the parameters. 

3.1. Eckart Potentials. The Eckart potentials are functions of the 
form 

ue x ve x 
= (1-e*) 2 ' 

where x is the physical distance variable, and u and v are parameters. 
Rewriting the Hamiltonian, 7i, in terms of z = e x we get 

HZ vz 

U = z 2 d zz + (1 - v)zd z + + -, (3.16) 

I — z {1 — z) A 

where the spectral parameter v has been encoded into the first-order 
part of the operator — see the preceding Section for a discussion of 
this matter. 

Let us now look for the connection between Eckart potentials and 
hypergeometric functions. The hypergeometric series, 

F(a,/3, T ,z) = F + F lZ + F 2 z 2 + ... , 

can be defined as the solution of the hypergeometric differential equa- 
tion: 



Q(F) = 0, F(0) = 1, 



(3.17) 
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where Q is the familiar second order differential operator 

Q = z{\- z)d zz + (7 - (a + /3 + l)z) d z - a/3. (3.18) 

In order to work with spectral residues we need to somehow convert 
the hypergeometric operator into an operator of the form z 2 d zz + lower 
order terms. An obvious way to do this is to left-multiply by zj (1 — z). 
To wit, set 

Qx 



l-z 



= z 2 d zz + (1 - u)zd t - ^— [(a + /3 + v)zd z + a(3] 

where v = 1 — 7 is the required spectral parameter. The upshot of all 
this is that one can equally well specify the hypergeometric series F as 
the solution of 

G 1 {F) = 0, F(0) = 1. 



In Proposition of the preceding section we noted that every 
second-order differential operator is gauge-equivalent to a Schrodinger 
operator. Following the method described there, we conjugate Q\ by a 
multiplication operator fi = e a where 



a{z) = - 
1 



(a + P + v) log(l - z) 



2 

li{z) = {l-zf a+(3+ ^l 2 . 

The end-result of the gauge-transformation is the Eckart Hamiltonian 
shown in ( p.!6| ). Using ( |2.12| ), an elementary calculation will show that 
the potential parameters are related to the hypergeometric parameters 
by 

u = -(a- p-v){a-p + v), (3.19) 
1 

v = -(a + P + v)(2-a-P — v). 

Turning next to spectral residues, we expand the Eckart potential in 
a series with respect to z: 



U = ^{u + kv)z k . 



In Proposition |2]^ we proved that the spectral residues of a Schrodinger 
equation are given by formula ( p. . 1| ) . From this formula and from the 
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above expansion we obtain the expression for the spectral residues of 
the Eckart potentials: 

p„{u,v) =n , (3.20) 

\p\=n 

where as usual the index ranges over all compositions p = (pi, . . . ,pi) 
of n. We have arrived at the following result. 

Proposition 3.1. In the even case, n = 2m, one has 

m—1 

p n (u, v) = K Y[ [{u + v + a n>i ) (u + v + a n>i+1 ) + b^v] , (3.21) 

i=0 

where K = [)p a n,i = i{n — i) and b ni = n — 1 — 2i. In the odd 
case, n = 2m + 1, one has 

m—1 

p n (u, v) = K(u + v + a n>m ) Yl [(« + v + a n>i ) (u + v + a n>i+1 ) + b^v] . 

(3.22) 

Proof. We saw above that the eigenfunctions of the Eckart Hamilton- 
ian are the product of a power of (1 — z) and of a hypergeometric 
series. Hence by the gauge-invariance of spectral residues as asserted 
in Proposition |2^ above, the spectral residues of the Eckart Hamilton- 
ian can just as well be calculated by taking residues of the coefficients 
of the hypergeometric series. 

These are of course the well-known: 

(3.23) 



" n\ [ T ]« ' 

where [x] n denotes the rising factorial x(x+l) . . . (x+n—1). Calculating 
the residue of F n at v = n (recall that 7 = 1 — v), and using the gauge- 
invariance of p n we obtain 



Pn{U, V) 

From ([3191) it follows that 



a nl(n — 1)! 

(a + i)(p + i)(a + n - 1 - i)((3 + n - 1 - i) = 



(3.24) 



[u + v + a n:i )(u + v + a nji+ i) + b^v , (3.25) 

' v=n 

Another calculation will show that if n is odd, say n = 2m + 1, then 

-(a + m)((3 + m) = u + v + a n>r7 
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This observation explains why the (— l) n factor in (|3.24 ) can be dis- 



carded. The proposition now follows immediately. □ 
To conclude the discussion of the Eckart potentials let us mention 



some identities that come about as a consequence of Proposition |3TT 



The first of these was mentioned in Theorem [L2| of the introduction. 
For the proof, we simply set u = in ( p.21|) and (|3.22| ) above, and note 
that 

{v + a n;i )(v + a ri:i+l ) + b^v = (v + i(i + l))(v + (n-i- l)(n - i)) . 

A similar factorization identity can also be obtained by setting u = —v 
in ( 3.22 ) and (13.21 ). In order to describe this identity, for positive 



integers k < n let us set 



|p|=n r 

e( P )=k 



Corollary 3.2. In the even case, n = 2m, the above numbers c n ^ are 
given by the coefficient of x k in 

m—l 

Y[ [in -l- 2ifx + z(z + \){n - 1 - i)(n - z')Y 

An analogous statement holds in the odd case, n = 2m + I, but with an 
additional factor of mim + 1) appended to the above expression. 

3.2. Poschll- Teller Potentials. Turning next to the Poschl- Teller 

potentials, these are functions of the form 

u v 
+ 



cosh 2 (x/2) sinh 2 (x/2)' 

where again x is the physical distance variable, and u, v are parameters. 
Rewriting the Hamiltonian, 7i, in terms of z = e x we get 

H = z% z + (1 - u)zd z + + ^-p. (3.27) 

Let us next consider how 7i can be related to the hypergeometric op- 
erator Q. To begin, let us rewrite the latter using ( as the independent 
variable: 

g = c(i - 0% + (7 -(« + /? + 1)0 <9c - <*P- 

We note that a change of variables cosh x = 2( — 1 will transform the 
leading term of Q into — d xx . Since we would like the leading term to 
be — z 2 d zz (the minus sign doesn't really matter), we employ z = e x 
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as the independent variable. To make the change of variables we note 
that 

(z + l) 2 Az 2 
li - ' ^~~ 2 

In this way we obtain 



C=^^ d c = -^-d z . (3.28) 



-G = z d zz + (1 - a - (3)zd z + 



(2 7 - 1)— — + (2a + 2(3 + 1 - 2 7 ) z<9 2 + a/5 

z+l z — 1 y 

The roots of the indicial equation are a and /3, and therefore, in or- 
der to obtain solutions that are monic power series we make a gauge- 
transformation: 

Gi = -z- a oQoz a 

= z 2 d zz + (1 + a - f3)zd z + 

+ ((2 7 - l)^y + (2a + 2(3 + 1 - 27)^-) (^ + a) 

Employing the abbreviations 

A = 7 - 1/2, (3.29) 

/x = a + /3+l/2-7, 

V — (3 — CK, 

Az /i^ 

turn = H - 

v ' z + l z-1 

the above operator can be more compactly rewritten as 

Gi = z 2 d zz + (1 - v)zd z + 2w(z)zd z + (A + \i - u)w(z) 

From there, the following gauge transformation will give the Poschl- 
Teller Hamiltonian: 

n = ((1 + z) x (i - zY) o g 1 o ((1 + z )-\i - z)~») , 

provided the potential parameters are related to the hypergeometric 
parameters by 

u = A 2 — A, v = fx — jj 2 . 

Turning next to the eigenfunctions of Ti we encounter a difficulty: 
we do not get a power series if we substitute ( = (z + I) 2 / (Az) into the 
usual hypergeometric series F(a, (3, 7; (). Note though that Az/ (1 + z) 2 
is analytic in z, and so we can overcome this difficulty by employing 
a solution of the hypergeometric equation that is analytic, or nearly 
so, in Just such a solution can be found on Kummer's list of 24 
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solutions to the hypergeometric equation — see Section 2.9 of 0. The 
solution to G{<p) = we require is 

(p(() = C a F(a, a + 1 - T , a + 1 - /3; C" 1 ). 
Let us summarize this observation in the following manner. 
Lemma 3.3. The unique series solution <fi = <p(z) to 

= 0, 0(0) = 1, 

is given by 



l + z) 



We are now in a position to factorize the spectral residues of the 
Poschl- Teller potential. Expanding the Poschl- Teller potential in a 
power series we obtain 



U{z) = y2 ku z k + ^2 kui 



k even k odd 

th 



where uo = v — u and u\ = v + u. Hence, by Proposition ^]5| the ri 
spectral residue, is given by 

Pn = nJ2^< P) < P \ (3-30) 

\p\=n 

where o(p) is the number of odd elements of a composition p = (pi, . . . ,pi), 
and e(p) is the number of even elements. 

Proposition 3.4. In the even case, n = 2m, one has 

p n = KY[(uj + 2k 2 uo + k 4 - k 2 ^j, (3.31) 

0<fc<n 
k odd 

where K = ,, 1 . In the odd case, n = 2m + 1, one has 

n\(n— 1)! ' ' 

p n = K Ul \\ (u\ + 2k 2 u Q + k A - k 2 ^j , (3.32) 

0<k<n 
k even 

Proof. Let <f>(z) = 1 + (fiiz + faz 2 + . . . be the unique series solution to 

0i(0) = O, 0(0) = 1. 



By Lemma 3.3 , if (fi + A — v)/2 = —j, j G N, then <f>(z) is (1 + 
z y-\-n times a jth degree polynomial of z/(z + l) 2 whose coefficients 
are polynomials of A, /i, v over (\ — v){2 — v) . . . (j — v). Hence for all 
n > j the residue of <p n at v = n must have \x + A + 2j — n as a factor. 
By the gauge invariance of spectral residues, the same must be true for 
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p n . The same reasoning also applies to show that A — fi — 1 — 2j + n 
is a factor of p n for all j — 0, 1, . . . , n — 1. 

We have now shown that function of A and fi, is some factor 

K times 

n-l 

+ A + 2j - n)(A - fi - 1 - 2j + n). 

3=0 

Now 

u = v — u = \i + A — fi 2 — A 2 , 
Mi = v + u = /i — A — fi 2 + A 2 = (A — + A + 1), 
and hence for all 

m 2 + 2k 2 u + k 4 -k 2 = 

(A + \i + k - 1)(A + fi - k - 1)(A - fi + k)(X - // - k) 

Using (|3.30|) , a degree count, and an examination of the coefficient of 
we see that K = n \^-iy. ■ The proposition now follows immediately. □ 

Here are some interesting identities arising from the preceding fac- 
torization. 

Corollary 3.5. The coefficient of x k in 

(x + 2 4 )(x + 4 4 )(x + 6 4 ) . . . (x + (2m) 4 ) 
is given by the following composition sum: 

^ n\ n! Yl^i 



p|=2m+l V V 
o(p)=2k+l 



(3.33) 



Proof. Upon setting x = u 2 , u = 1/2, the above identity follows from 
the odd, n = 2m + 1, case of Proposition [T4L □ 



If we use the same summand as above, but restrict ourselves to com- 
positions without even numbers we obtain the coefficients of a closely 
related product. 

Corollary 3.6. The coefficient of x k in 

(x + 2 4 - 2 2 )(x + 4 4 - 4 2 )(x + 6 4 - 6 2 ) . . . (x + (2m) 4 - (2m) 2 ) 
is given by the following compositions sum: 

n\ n\ 

<5n St\' 



Enl nl T-r 
ll Pi 
Q__ <J_, J. J. 



p|=2m + l P P 
o(p)=2fe + l 
e(p)=0 
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Proof. Upon setting x = uf, Uq = 0, the above identity follows from 



the odd, n = 2m + 1, case of Proposition 3.4. □ 



Identities featuring k 4 and k 4 — k 2 , with k odd, can be obtained in an 
analogous manner. 

3.3. Morse Potentials. Finally, let us consider Morse potentials; these 
are functions of the form 

ue x + ve 2x , 

where x is the physical distance variable, and u, v are parameters. 
Rewriting the Hamiltonian, 7i, in terms of z = e x we get 

H = z 2 d zz + (1 - v)zd z +uz + vz 2 , (3.34) 

The eigenfunctions of the above operator can be related to confluent 
hypergeometric functions by a gauge transformation. The confluent 
hypergeometric series can be defined as the solution of the confluent 
hypergeometric differential equation: 

j{4>) = o, 0(0) = l, 

where J is the second order differential operator 

J = zd zz + (7 - z)d z + a 

We need to introduce an extra scaling parameter, u, stemming from 
the change of variables z 1—* — 2uoz. Now, the relevant operator is given 
by 

J = zd zz + (7 + 2uz)d z + 2ua, 

and the solution, 

(f>(a, 7, LO] z) = 0o + 4>xz + 4>2Z 2 + ... , 

of J[(f\ — depends on 3 parameters. 

The equation J7"[0] = is equivalent to a certain relation between <fr n 
and n +i, and yields the following modification of the familiar formula 
for the coefficients of 0: 

(-2) n Lu n [a} n 
0n ~ n\ [7]" • 

Proceeding as we did in the case of the Eckart potentials we modify 
J so as to obtain an operator with leading term z 2 d zz : 

Ji = zj = z 2 d zz + (1 — v)zd z + 2uoz 2 d z + 2toaz, 
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where again we employ v — 1 — 7. Next we convert J\ to (nearly) self- 
adjoint form by means of a gauge transformation with gauge factor 
p = e a , where 



(J — I UJ = UJZ. 



Hence, p = e ulz . The result is the Hamiltonian operator shown in ( |3.34| ). 
where the potential parameters are related to the hypergeometric pa- 
rameters by 

u = u(y + 2a-l), (3.35) 
v = —uj 2 . 

To obtain the corresponding spectral residues we again employ (|1.1| ): 
p n (u,v) = n ^2 



u °(p) v e (p) 



S p S p i 

\p\=n 
Pi G{l,2} 

where the index now ranges over compositions p = (pi, . . . ,pi) of n 
consisting of 1 and 2 only, and where it is necessary to recall that o(p) 
denotes the number of odd elements of a composition p = (p 1; . . . ,p t ), 
and e(p) the number of even elements. 

Proposition 3.7. In the even case, n = 2m, one has 

p n (u,v) = Kl[(u 2 + k 2 v), (3.36) 

0<k<n 
k odd 

where K = ., 1 . In the odd case, n = 2m + 1, one has 

n\(n— 1)! ' ' 



Pn[U : V) 



KuY[(u 2 + k 2 v), (3.37) 



0<fe<n 
k even 



Proof. We proceed in the same way as in the proof of Proposition 
(|3"TT|). By the gauge- invariance of spectral residues, p n is equal to the 
n th spectral residue of the confluent hypergeometric operator. Taking 
the residue of n at v = n yields 

2 n uj n [a] n 
(n\) 2 ' 



Next, using Q3.35D we see that in the first case, where n = 2m, k = 2i+l, 
we have 



u + k 2 v 



Aoj 2 (ol + m + i)(a + m — i — 1). 



SPECTRAL RESIDUES 22 

In the second case, where n = 2m + 1, k = 2i we have 
u 2 + k 2 v 



u 



Auj 2 (a + m + i)(a + m — i), 
2cu(a + m). 



The proposition now follows immediately. □ 

Setting v — 1 in the above Proposition immediately yields the following 
identities: 

Corollary 3.8. The coefficient of x k in 

(x + 2 2 )(x + A 2 )(x + 6 2 ) . . . (x + (2m) 2 ) 
is given by the following composition sum: 

y-v n\ n\ 

|p|=2m+l 1 v 
o(p)=2fc+l 
Pi £{l,2} 

T/ie coefficient of x k in 

(x + l 2 )(x + 3 2 )(x + 5 2 )...(x+ (2m - l) 2 ) 
is given by the following composition sum: 



|p | — Ira 
o(p)^2fc 

Pie{i,2} 



5p ^p' 



4. The Darboux transformation 

We begin this section by describing the Darboux transformation 
for linear, second-order differential equations. Let us consider a self- 
adjoint, second-order operator 

H = d xx + U(x), 

and ask whether it can be factored into a product of first order op- 
erators. It is not hard to see that the answer is "yes", and that the 
factorization must be of the form 

n = (d x -w)(d x + w), 

where the function W(x) — we will refer to it as a prepotential — is 
related to the potential by U — W x — W 2 . Reversing the order of the 
above factorization we obtain another self-adjoint operator: 

H = (d x + W)(d x -W) = d xx + U, 
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where now U = —W x — W 2 . Knowing U one has to solve a Ricatti 
equation in order to obtain W. Alternatively, one can describe the 
prepotential as W = —d x (logip), where ip(x) is a solution of 7i(ip) = 0. 

The process of going from U(x) to U (x) is called the Darboux trans- 
formation. It is evident, that at the level of the prepotential, the Dar- 
boux transformation corresponds to a negation, and that therefore the 
Darboux transform of U takes us back to U. 

There is a close relationship between the (formal) eigenfunctions of 
TC and H. Indeed, one can write 

(d x + W)(d x - W)(d x + W) = (d x + W)H (4.38) 

= H(d x + W), 

From here, suppose that if) is a formal eigenfunction of H, i.e. that 
= -V0. Setting ip = ip x + Wip we employ (|4. 38|) to conclude that 
H$) = Xi>. 

In order to apply the Darboux transformation in a truly analytic 
setting, one has to take into account factors like square-summability 
of the eigenfunctions, and smoothness of the potential. Here, we are 
more interested in formal properties of the Darboux transformation. 
Moreover, our approach is based on power series in z = e x , rather than 
on functions of x. 

Let us therefore adapt the notion of a Darboux transformation to 
the present setting. We consider a second order operator 

n = z 2 d zz + (i - u)zd z + u, 

where U(z) is a formal power series with vanishing constant term. We 
do not assume that the coefficients U n commute. Owing to the fact 
that U(z) has zero constant term, H cannot be factored into a product 
of first order operators. In order to obtain a factorization we must add 
a constant. To wit, 

U + zy2/4 = (zd z - u/2 - W)(zd z - v/2 + W), 

where W(z) = W\z + W 2 z 2 + . . . is another formal power series related 
to the potential by 

U = zW z - W 2 . 

At the level of coefficients this is equivalent to an infinite number of 
graded relations: 

n-l 
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Note that W n = n _1 U n + lower order terms, and hence the above 
relations can be recursively inverted, so that if we know the coefficients 
of U, we can calculate the coefficients of W. 

Reversing the above factorization we arrive at the partner operator 

U = z 2 d zz + (1 - v)zd z + U, 

were U = —zW z — W 2 . Here are the first few relations that describe 
the formal Darboux transformation: 

Ui = -U h U 2 = -2U?-U 2 , U 3 = -2U 1 3 -2U l U 2 -U 3 

It does not seem possible to obtain a closed sum formula for the above 
relations. 

Next, let <j){z) = 1 + <p\z + <p 2 z 2 + . . . be the formal monic solution 

of 

H{4>) = 0. 

We saw in Section ^| that the coefficients of will be certain rational 
functions in v. We set 

4> = z<p z + (-z//2 + W)(f>, 

and note that since 

(H + u 2 /A)(zd z - u/2 + W) = (zd z - u/2 + W)(H + v 2 /A), 
we can conclude that 

H{4>) = 0. 

We are now ready to prove Theorem |1.3| . Note that 

ij) = -2v- l ~$ 
is the unique monic solution of 

W(V0 = o. 

Hence the n th spectral residue of U, call it p n is given by 



Since 

it follows that 



2 

p n = Res(0 n ; u = n). 

n 



(n — v/2)cf) n + lower order terms, 



2 / n 



-n\ n -2) pr ' 



where p n is the n th spectral residue of U. The theorem is proved. 

Finally, let us see how the relation between the spectral residues 
and the prepotential leads us to the inversion formulas ( |1.6| ) and ( |1.7|) . 
With Theorem |0| at our disposal it is not difficult to describe W in 
terms of the p n . 
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Proposition 4.1. 

Wn |^ (Pi + V2){.V2 + Pz) ■ ■ ■ (Pl-l + Pi) 

t(p) odd 

Proof. Note that 2zW z = U - U, or equivalently that W„ = (£/„ - 
f/ n / (2n). Furthermore, Theorems |1.1| and |1.3| imply that 

Un= y „ 

ii (Pi + P2XP2 + Pa) • • • (pz-i +Pi) P 

|p|=n 

(Pi + P2XP2 + Ps) ■ ■ ■ (pi-i + Pi) Pp 

\p\=n 

The desired conclusion now follows immediately. □ 

The inverse relation, i.e. a formula for the spectral residues in terms 
of the W n is shown in equation (|1.7[). Let us prove it now. Let and 
denote, respectively the monic solutions of 7^(0) = and TC((p) = 0. 
Alternatively, we could consider and together as the solutions of 
the following coupled first-order equations: 

z<t>z-^<P + W<l> = ~ij> 

0(0) = 0(0) = 1 
Introducing sum and difference variables: 

<x = (0 + 0)/2, 5 = (0-0)/2, 

we obtain 

ZC r z = -W5, cr(0) = 1, 

z8 g -u8 = -Wo, 8(0) = 

In terms of power series coefficients we have 

na n = -W n -i8\ - ... - Wi8 n -! 

(n - V)8 n = -W n - Wn-lCl - • • • - Wi(T n -i 

Now 

Res(0 n ; v = n) = p n , 



and by Theorem |1.3| , 

Res(0 n ; v = n) = -p n . 
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Consequently, 

p n = Res(5 n ; v = n). 

It is evident that given W\, W2, ■ ■ ■ the above relations can be solved 
for a n and It is also evident that each <j n will be a sum of products of 
even numbers of W n and that each 8 n will be a sum of products of odd 
numbers of W n . Indeed, solving the above relations for 5 n we obtain 

Sn= E A x .„ x Jt K ik x ik x 



p| — 71 

e(p) odd 



v — n (— s 4 ) v - s 3 (— s 2 ) v — s\ 



where as usual p = (pi, . . . ,pi) is a composition of n, and = pi + 
. . . + pj. Taking the residue of the right hand side at v = n, equation 
Ql.Tf) follows immediately. 

5. Applications and Conclusion 

The above discussion indicates how the notion of a spectral residue 
can be used to give a formal series treatment of differential equations 
with spectral parameter. Our focus has been on second-order, one- 
dimensional problems. This is a well-understood theory with a large 
selection of techniques for exact solutions. It has been our aim to 
demonstrate that such techniques lead to useful information once we 
make the jump to the formal domain. The author feels that the ap- 
plication of the spectral residue methodology to more general classes 
of eigenvalue problems, e.g. non-linear and higher-dimensional linear 
problems, will yield additional combinatorial results. 

In conclusion, let us also indicate some applications for the identities 
and the inversion formulas developed here. We will consider two topics: 
enumeration and integrable systems. An application to shape-invariant 
potentials has been reported on in [15]. The author hopes that the 



existence of such applications will serve as an incentive for a further 
investigation of spectral residues. 

5.1. Enumeration problems. Consider a selection of k triples of 
numbers (a, b, c) with l<a<b<c<n. Alternatively we could talk 
about a 3 x k rectangular tableaux whose entries are numbers from 1 
to n (repetitions permitted) with the constraint that the columns be 
arranged in ascending order. A natural question in this setting is the 
following. If the k triples are selected in a uniformly random manner, 
what is more likely: a multiple occurrence of a number in the second 
row, or a multiple occurrence in the third row. The answer is that 
duplication is more probable in the third row (or indeed in the first 
row by consideration of symmetry) than in the second row. 
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This fact is straightforward consequence of the identity shown in 
1.4). To see why let us remark that ( m 2 _1 ) is the number of triples 



with m as the third entry. Therefore the generating function that 
counts tableaux with distinct entries in the third row is 

. . / 1-2 \ / 2-3 \ / n(n-l) 
F(x) = 1 + — - x 1 + — - x ... 1 + — '- x 



This is more or less the right hand side of equation ( p..4|) . Next we 



note that m(n — m) is the number of triples with m as the middle 
entry. Therefore the generating function for the number of tableaux 
with distinct entries in the middle row is given by 

n-1 | | 

Y[ (1 + m{n - m)x) = ^ — — x"^ (p) 

m=l | p |=„ Sp Sp ' 



This is closely related to the left hand side of ( |1.4| ). Indeed, rewriting 
(11.41) we have 



F(x) 




n-e(p) 



The weights YliPi/(2 Pi ~ 1 ) m t ne above right-hand side are less than 
one, and the desired conclusion follows. 

The other identities derived in Section |3| have similar enumerative 
interpretations. The gist of these interpretations is that various mea- 
sures that count sparseness of triples relative to their largest entries 
correspond to certain other measures that count sparseness relative to 
the middle entries. These developments will be reported on elsewhere. 

5.2. Integrable systems. It has already been stated that spectral 
residues are meant to furnish a kind of spectral information for formal 
differential equations with a parameter. It is natural then to enquire for 
a physical interpretation of the spectral residue sequence. We will show 
that it is makes sense to regard the spectral residues as a formalized 
version of the potential's scattering data. 

Such an interpretation comes about when we inquire about the evo- 
lution of the spectral residues concomitant to the evolution of the po- 
tential function with respect to the KdV equation. 

Proposition 5.1. Suppose that the potential series U(z,t) evolves ac- 



cording to the non- commutative version of the KdV equation (see [16| 
for a list of references): 

U t = D 3 (U) - 3D(U)U - 3 UD(U), 
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where D = zd z . Then, the spectral residues p n evolve according to the 
rule 

d t (pn) = n 3 p n . 

The proof is straightforward; it relies on the following identity 

(pi + ...+Pnf ~ {pf + ...+Pn) = 

53 (Pi + • • ■ +Pj)(p j+ 1 + . . .+Pn)(Pj +Pj+l)- (5.39) 

5 

Thus if we interpret the n in p n as square roots of the energy ( see the 
discussion following equation ( |2.11[) ) we see that the spectral residues 



evolve according to the same rule as the scattering data (reflection and 
normalization coefficients) of KdV wave profiles. 

We now have the following quite natural interpretation of the trans- 
formations relating the coefficients of U(z) and the spectral residues, 
i.e. equations (1.1) and (1.2): theses inversion formulas describes a 



formalized version of the well-known inverse scattering transformation 



for the KdV equation [[H]]. As such, the soliton solutions correspond 
to finitely supported residue sequences, and indeed it is possible to re- 
cover the well known n-soliton formula [|18|, [19] by means of equation 

(O- 

It must be noted that the formula (11.2) and its relation to solitons has 



been established previously by Wadati and Sawada [HJ. The identity 
Q5.39Q also occurs in that paper. Indeed, these authors went so far as 
to derive the Gelfand-Levitan-Marchenko equation from ( |1.2| ). This 
is not surprising, in light of the above interpretation of (|1.2|) as the 
restorative half of the inverse scattering transform. The missing idea 
in Wadati and Suwada's article was equation (|1.1| ): the forward portion 
of the formalized I.S.T. 

The next logical step in this development is the representation of 
other integrable systems by means of combinatorial inversions. It is 
reasonable to expect that the spectral residue formalism will be useful 
in translating the spectral problems that represent integrable equations 
into appropriate analogues of ( |1.1|) . The possible benefits of such a pro- 
gram are two- fold. First, the translation of scattering transforms for 
other integrable systems may yield new combinatorial results on par 
with the results for the Schrodinger equation reported on in this arti- 
cle. Second, by analogy with harmonic analysis it may be feasible to 
characterize classes of solutions to KdV by means of spectral residues. 
We already know that solitons correspond to finitely supported residue 
sequences. Is it possible to characterize compactly supported or expo- 
nentially decaying solutions in terms of an appropriately rapid decay 



SPECTRAL RESIDUES 



2!) 



of the residue sequence? If such an approach is successful, it seems 
reasonable that it could be carried over to other integrable equations. 
These matters call for further investigation. 
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